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In this paper, we report on the sharp peak observed in the third harmonic voltage response generated by a bias 
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1. Introduction 

Harmonic signals created by ac excitation have 
been used for a long time in the investigation 
of high temperature superconductors (HTSC) |Q. 
Particularly, a third harmonic signal is a conse- 
quence of non-linear effects and its measurement 
is a powerful technique to investigate, and more- 
over, to discriminate between the different pos- 
sible origins of these non-linearities. Different 
techniques of third harmonic measurement exist 
for investigating properties of HTSC devices espe- 
cially in microwave range 0-|6), or for character- 
izing material properties such as thermal conduc- 
tivity [^j-pUl, specific heat pL current response 
|l2[ , or critical temperature |13[]. 

The existence of a sharp maximum in the third 
harmonic can be easily understood from a qual- 
itative point of view. At low temperatures the 
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voltage response of the superconducting phase is 
almost zero. On the other hand, in the normal 
phase we have with high accuracy an ohmic be- 
haviour and the response at higher frequency har- 
monics is again small. A sharp maximum related 
to a viable variety of nonlinear phenomena is ex- 
pected only close to the critical temperature. 

As an example, measurement of the third har- 
monic of the ac susceptibility revealed the role 
of interplanar spacing on formation of flux lines 
and demonstrated that the irreversibility field in 
Bi2212 is determined by surface barriers [Jl3| , jl4| . 
This is related to vortex motion, which is one 
of the possible origins of the non-linear effects 
in electrical properties in HTSC See also 

references |l^] and flr7[ where a mechanism is 
discussed for the generation of third harmonics 
based on the anisotropic scattering rate of normal 
charge carriers included in the Boltzmann equa- 
tion. 

One of the other possible roots of this non- 
linear behaviour is given in ref. Jl8[ , where a 
theoretical derivation is presented to explain non- 
linear effects in electrical conductivity close to the 
transition temperature. From this paper ]l8| ], it 
is expected that the measurement of the third 
harmonic voltage allows getting hold of supercon- 
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ductor constants of the BCS theory. In addition, 
the authors Q , []l9) , investigate theoretically the 
third harmonic voltage generation produced by 
the thermal response of the film to a harmonic 
current. Third harmonic response is definitely 
related to different roots of non-linear effects and 
should allow one to gather new information about 
the mechanism of HTSC non-linearity. 

In the present paper, systematic measurements 
of the voltage response to a purely harmonic 
current crossing HTSC strips patterned in a 
YBa 2 Cu307-,5 (YBCO) thin film and in two 
Las-zSr^CuCU (LSCO) thin films with x=0.15 
were performed as a function of temperature. 
First and third harmonics were recorded for the 
whole temperature range, from room temperature 
to a few Kelvin below the normal-superconductor 
transition. A sharp peak was observed in the 
third harmonic voltage close to the critical tem- 
perature, exhibiting a lambda-shaped (A-shaped) 
behaviour. This A-shaped peak was investigated 
for different current amplitudes and frequencies. 

Sample and measurement descriptions are 
given in part 2 of this paper. Complete experi- 
mental data processing and observations are pre- 
sented in section 3. In section 4, mechanisms of 
third harmonic signal generation are briefly dis- 
cussed. 

2. Experimental setup 

Measurements were first performed on two dif- 
ferent strips designed in a YBCO thin film, 180 
nm thick, deposited by DC sputtering on a MgO 
substrate, as described in [^0| and pi] . The 
smallest is 1.59 mm long and 77 fim wide, and 
the largest is 2.87 mm long and 474 /im wide. 
Both were chemically etched in acid solution. A 
second set of measurements was performed on 1 
mm long and 100 /xm wide strips patterned on 
two different LSCO films, LSCOl and LSC02, 
both deposited by sputtering on a SrTiOs sub- 
strate. The first one was synthesized at 850°C, 
under 1.83 mbar of pure oxygen, and annealed 
under pure oxygen at 10 mbar at 850°C for half 
an hour. For the second one, the temperature of 
deposition and annealing was decreased to 840° C. 

The samples were placed in a cryostat and 




Figure 1. Scheme of the electrical measurements. 
The current in the superconductor sample (sym- 
bolized by its temperature-dependent resistance 
R(T)) is governed by the voltage generator and 
the load resistance Rl- The voltmeter across 
Rl allows to monitor the current in the sample. 
The rms voltage harmonic response (V rlLU -r */2, 
n = 1, 3) across the sample is measured by a 
lock-in amplifier and recorded as a function of 
temperature. 



exposed to a helium flux. The connection be- 
tween the electronic equipment and the strip was 
made via coaxial cables in a four-points geome- 
try (fig. |l|). Strips were bias current, the latter 
being sinusoidally modulated with amplitude Iq 
and frequency / = uj/2tt. The rms current value 
Iq/ -\/2 was monitored during the whole experi- 
ment by measuring the voltage across the load 
resistance Rl, chosen large enough to have negli- 
gible variation of the current (less than 10 %). Us- 
ing a digital lock-in amplifier (model SR830 from 
Stanford Research Systems), first (Vi u ) and third 
(y^) harmonics of the voltage developed in the 
sample were successively recorded simultaneously 
with temperature and time. 

As was mentioned in previous papers 0, Q, 
the current has to be much smaller than a maxi- 
mal value Iq : 

where R' is the first derivative of the electrical 
resistance with respect to temperature and G is 
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the thermal boundary conductance between the 
strip and the substrate, which is maintained at 
the temperature of the helium bath. G can be 
expressed through the thermal boundary conduc- 
tivity g, the length L and the width w of the 
strip: 



G = gLw 



(2) 



In the case of the YBCO film, calculating Iq 
from standard values of g (2000 W/Kcm 2 ) and 
R' (1 O/K) in the normal state, we found the 
maximal current to equal 1.6 A for the narrowest 
strip and 5.2 A for the widest one. In our case, 
the maximum current used had an amplitude of 
7 niA, far enough from the limit given by Eq. (Q). 
Concerning the LSCO films, the biggest current 
amplitude used was 1.4 mA. 

3. Experimental observations 

3.1. Experimental Data processing 

Before showing our results, we explain in some 
detail the experimental data processing, taking 
as an example the measurements obtained from 
the YBCO film. The goal of this processing is to 
calculate the first and the second derivatives of 
the resistance with respect to temperature, which 
we need to estimate the thermal conductance G, 
as described in ref. p9j. Since our cryogenic 
system induces strong temperature noise, tem- 
perature data were smoothed over time, using a 
Savitzky-Golay filter j2^| , which performs a local 
polynomial regression to determine the smoothed 
value for each data point. In our case, 2 nd order 
polynomial regression was used, taking into ac- 
count 25 points to the left and 25 points to the 
right. 

The same smoothing was used to filter the elec- 
trical resistance of the strip, R. Since the fre- 
quency lu is in the kHz range, we consider the 
superconducting strip as a purely ohmic device 
where inductive and capacitive components are 
negligible. The resistance is then deduced from 
Viu by Ohm's law: 

fi> - Vl " 
Jo 

In this case, only 5 points to the left and 5 points 
to the right were taken into account to perform 
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the local regression, to avoid any big distortion 
close to the superconducting transition. Exactly 
the same procedure was applied for smoothing the 
third harmonic voltage. 

Thus, sets of time-dependent smoothed-data 
were obtained and used to calculate the first and 
second derivatives of the temperature (T, T) and 
resistance (R, R) with respect to time t, differen- 
tiating using Savitzky-Golay smoothing. For this 
step, the local smoothing was performed on the 
smoothed-data sets. The local regression taking 
into account 21 points for the temperature set and 
11 points for the resistance set, allowed to deduce 
the first and the second derivatives for each data 
point. Thus, it was straightforward to determine 
for each point in time, the values of the first and 
second derivatives of the resistance with respect 
to temperature T: 



R'(T) 



R 
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RT-RT 
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All the steps of the experimental data processing 
are gathered in figures || and [|, presenting, respec- 
tively, the temperature, the resistance and their 
respective derivatives obtained for the narrowest 
strip excited by a current of 4 mA rms. The tem- 
perature was shifted by a few Kelvin during the 
measurements due to the poor thermal conductiv- 
ity of the sample holder on which our temperature 
diode and sample were glued; we then shifted all 
temperature data sets by 4 K in order to have a 
reasonable T c of 90 K for our YBCO sample. 

Figures [| and || give the resistance and its first 
and second derivatives with respect to smoothed 
temperature, deduced from the time-dependent 
data sets shown in figures || and ||, with use of 
eq. |I| These results were used to calculate the 
thermal conductance G between the strip and the 
heat sink, as described in jl9| and to evaluate the 
contribution of the thermal effect in the peak ob- 
served in the third harmonic voltage Vs^,, which 
will be described in the next section. 

3.2. Lambda shape of the 3 rd harmonic 
voltage V^u 

Figure || gives the temperature dependence of 
the third harmonic response of the narrowest strip 
patterned on the YBCO film to the current exci- 
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tation with amplitude Iq equal to 5.65 mA and a 
modulation frequency of 1024 Hz. 

To have an appreciation of the common be- 
haviour of the third harmonic response for dif- 
ferent samples, we performed the same kind of 
measurement on two LSCO films. Figures [?] 
and H present for each tested sample (LSCOl 
and LSC02) both the resistance and the third 
harmonic response obtained for currents ranging 
from 10 fi A rms to about 1 mA rms, as a func- 
tion of the temperature. The sharp transition of 
the resistance presented in figure ^a gives rise to 
a unique lambda-shaped peak in figure ^b. In 
the case of a sample exhibiting a "double" transi- 
tion in the resistivity as for LSC02 (see fig. ||a), 
we observe two peaks in the third harmonic volt- 
age response, fig. |^b. This puts into evidence the 
high accuracy of this technique of measurement, 
firstly, to determine the critical temperature of a 
material in a very simple way, and further, to dif- 
ferentiate between different phases in a material 
and to give for each of them the critical temper- 
ature. 

The behaviour of the shape and the amplitude 
of the peak with respect to current characteristics 
has been explored. As we can see on figures [7)3 
and ||b the amplitude of the peak increase with 
the current. Figure [| presents the dependence of 
the peak value of the third harmonic response on 
the current amplitude in a log- log plot. A power 
law dependence of the third harmonic peak am- 
plitude on the current intensity is exhibited for 
each curve. But no clear power law emerges, the 
power ranging from 1.1 for the LSCOl film to 3/2 
for the YBCO film. We are also not aware of a 
theoretical prediction of a power law for the peak 
amplitude versus current. At T c , various contri- 
butions to the third harmonic play a role and are 
difficult to estimate. 

For T > T c , the situation is easier to analyse 
theoretically Jig ], p9| and a power law with an 
exponent 3 is predicted for the dependence of the 
third harmonic amplitude on the current. In our 
experiments, however, a complication arises due 
to the presence of a parasite current coming from 
our generator, which produces a third har- 
monic signal that plays an important role in this 
temperature range. This effect hides the real am- 



plitude of the third harmonic signal. Following 
the theoretical works @, |jl9), we have adopted 
a fitting procedure which takes into account this 
parasite current. In any case, even with these dif- 
ficulties, it is important to note that the third har- 
monic signal presented in this paper follows qual- 
itatively the expected lambda-shaped behaviour 
predicted in Jl9| . 

3.3. Experimental Data Analysis 

Figure [lO] presents the first attempt to con- 
nect theory and experiment. It shows the fit of 
the lambda shape, in the case of the 1 mA rms 
current excitation, using two different functions. 
The first one takes into account only the thermal 
effect. In the second one, both the thermal effect 
and the nonlinear current generated by the gen- 
erator, are considered. 

To have an appreciation of the reliability of the 
fit, the thermal resistivity between the film 
and the heat sink has been calculated from the 
thermal resistance Rh given by the fit, consid- 
ering the surface of the patterned strip (100 fim 
x 1 mm). If we consider the case where both 
the parasitic current 1^ and the thermal con- 
ductance effects are pesent, the fit gives a value 
of Rh equal to 580 K/W, which corresponds to a 
thermal resistivity equal to 5.8 x 10~ 5 m 2 K/W 
at the considered temperature (about 30 K). 
To our knowledge, no data has been published 
concerning the thermal boundary resistivity of 
LSCO thin films at this temperature. Then we 
compare the obtained value to the one theoret- 
ically predicted in |23| for a YBCO thin film 
deposited on MgO. The value predicted in the 
latter Q is about 10~ 6 m 2 K/W at 30 K. Our ex- 
perimental value is not far from this theoretically 
predicted value (about one order of magnitude 
bigger). Because the nature of substrate and film 
are different, some difference is to be expected. 
Furthermore, in the range of temperatures con- 
sidered (20 K < T < 60 K), the variation of 
the thermal boundary resistivity is rather large. 

The dependence of the peak on the current 
frequency was studied at fixed amplitude current 
(Iq = 1.41 mA). No differences in the peak were 
observed for frequencies ranging from 114 Hz to 
20.6 kHz. 
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4. Discussion and conclusion 

The main qualitative property of the third har- 
monic response is the sharp A-peak close to the 
critical temperature. Taking into account the 
third harmonic signal generated by the thermal 
oscillations in the film and by the parasite cur- 
rent from the generator, the experimental data 
have been fitted. From this procedure, we found 
an acceptable value of the thermal resistivity, , 
and a reasonable value of the parasite current, 
Thus, we have for the moment a qualitative 
agreement between experiments and theory. 

A more detailed investigation, beyond the par- 
tial analysis presented in section 3.3, of the de- 
pendence of the height and width of this peak 
on the current and geometry of the strip (espe- 
cially the thickness) will be the theme of further 
research. This should be helpful to discriminate 
the different contributions to the third harmonic 
response above T c , close to the critical temper- 
ature, i.e. the thermal oscillation part and the 
non-linear conductivity part. 

Other third harmonic sources (including the 
non-linear conductivity) would have to be taken 
into account to describe more accurately the 
whole peak. For example, a temperature shift 
of the peak is not yet explained by the proposed 
theories Q, Also, slightly above T c , we ob- 
serve a disagreement between our fit and the ex- 
perimental data in the tail of the peak. Let us dis- 
cuss in brief some mechanisms for the generation 
of harmonics in different temperature regimes and 
physical situations, which would have to be taken 
into account to describe the whole peak with a 
better accuracy. 

At very low temperature, in the vortex pin- 
ning regime or in a vortex-free situation, mod- 
ulation of the superconducting order parameter 
creates dissipation and harmonic generation [p4[ 
|25| ; the non-linearity is created by the current 
dependence of the superfluid density which can 
be easily described within the Ginzburg-Landau 
mode l p(| . Non-linear response of the vortex 
state |15|] can also create higher harmonics. How- 
ever, the complexity of the pinning mechanisms 



and the generation of a magnetic field by the ap- 
plied current makes it very difficult to develop a 
detailed theory for the low-temperature tail of the 
A-shaped maximum. For example, small super- 
conducting grains in a diamond anvil juj, 
or scanning of the surface of a thin film by a 
terminated coaxial cable [Q] require the solution 
of different electrodynamics problems. In our 
paper, we used a narrow strip and the vortices 
are created near the edges of the sample by the 
applied current. The generation of vortices can 
be significantly reduced if a Corbino geometry is 
used; it is necessary that a superconducting film 
terminates a coaxial cable. 

For temperatures slightly above the critical 
one, there is significant fluctuation conductivity. 
Usual ohmic conductivity of normal carriers is 
perfectly linear but fluctuation conductivity is 
not and can create an observable part of the har- 
monic generation |l8| . 

For higher temperatures, temperature oscilla- 
tions can create harmonics in the voltage response 
§], |§, HI, H). It is this mechanism which has 
been given most attention in our theoretical anal- 
ysis of the experimental data for T > T c (see Fig. 

For completeness, we would like to mention 
that analogously to temperature oscillations well 
above T c , a quite universal mechanism for 3 rd har- 
monic generation at microwave frequencies exists 
for the superconducting state below T c . The den- 
sity of superfluid particles depends on the applied 
current and finally the kinetic inductance of the 
superconducting sample depends on the current 
I = Iq cos(wt) 



L{I) 



Lo + LJ 2 . 



(5) 



Lq gives the linear contribution to the kinetic in- 
ductance and Li the non-linear one (see ref. pH , 
PI). Then, the general formula for the voltage 
response given by 

V(t)=Hl)f t (6) 

allows to calculate the third harmonic response 
El, §3: 



(7) 
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It appears in this formula that the third harmonic 
response increases linearly with the frequency. 
This explains that it is typically observed in the 
microwave range devices. In our range of frequen- 
cies (up to 20 kHz), this response is not expected 
to be observable. 

Apart from studying the interface heat re- 
sistance, for the fundamental physics of super- 
conductivity of layered cuprates it will be very 
interesting to investigate the bulk non-ohmic 
conductivity related to the kinetics of the normal 
phase. For further theoretical discussion see for 
example @, @ and [§§. 

In conclusion, the investigation of the A-peak of 
the third harmonic response is a good qualitative 
method for the detection of the superconducting 
phase transition. Because of the simplicity of the 
measurement, it can become a standard method 
in the materials science of superconductors and 
also in the industry for testing the quality of 
the thin films and cables. In order to separate 
the bulk and interface mechanisms for the third 
harmonic generation it is necessary to perform 
detailed investigations of superconducting films 
with different thicknesses. In this way applied 
research related to the development of supercon- 
ducting bolometers can become an indispensable 
tool for the investigation of the fundamental ef- 
fects in cuprate superconductors (29) and |30| . 
Our work demonstrates that the maximum of the 
A-shaped peak as a function of temperature can 
be very easily observed, even without any special 
care for signal cleaning. 

Because the voltage generator produces some 
third harmonic signal, if we wish to investigate 
the third harmonic far from T c , we have to com- 
pensate the parasite 3w-signal not only at the 
level of the data analysis, but also electronically. 
However, a very simple filter will lead to a re- 
duction of the ratio I^uj/h (~ 1/500) by a factor 
of about 10 only, so that the remaining parasite 
contribution must still be taken into account in 
the data analysis, as described in the theoretical 
works (equations (7.10) from Jl8| and (3.8) from 

Detailed investigation of the third harmonic 
response provides a reliable method for the in- 
vestigation of superconductor thermal properties 



like thermal conductivity, thermal boundary re- 
sistance and heat capacity. Moreover, if the 
thermal effects are carefully subtracted, the in- 
vestigation of the third harmonic signal will 
provide the opportunity to extract the relaxation 
time of the superconducting order parameter 
from the isothermal electric non-linearity. If a 
Corbino geometry does not create significant 
heating of the sample, this can be expected to be 
the best experimental setup for the observation 
of Cooper-pair depairing non-linearity. 
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Figure 2. Recording of temperature T vs time Figure 3. Recording of resistance R vs time (a), 

(a) , and first (b) and second (c) derivatives of its and first (b) and second (c) derivatives of its 

smoothed values T s . smoothed values R s . 
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Figure 4. The resistance and its first and second 
derivatives with respect to temperature deduced 
from eq. |I] and data sets extracted from figures ^ 
and|. 



Figure 5. The resistance and its first and sec- 
ond derivative with respect to temperature, on 
the temperature range indicated by the hatched 
areas on figures || and ||. 
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Figure 7. Measurements on LSCOl. a: Depen- 
dence of the strip resistance on temperature ; 
b: Dependence of the third harmonic signal on 
temperature for different current amplitudes. 



Figure 8. Measurements on LSC02. a: Depen- 
dence of the strip resistance on temperature ; 
b: Dependence of the third harmonic signal on 
temperature for different current amplitudes. 
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Figure 6. A-shaped temperature dependence of the third harmonic voltage Vs u for YBCO. 
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Figure 9. Dependence of the maximum of the third harmonic signal on current amplitudes for the three 
tested samples. 
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Figure 10. Fit of the lambda-shaped peak in the case of a 1 mA rms current measurement on LSCOl, 
taking into account the thermal effect only (with the thermal resistance Rh between film and substrate 
as parameter), or both the thermal and the parasite 1^ current effects. The equation used for the fitting 

j t r nr , RhRR'In , 5fl? R(2 R' 2 + RR" ) I* 

procedure is: V 3u = Rhu> + - J ^— A — ~ H — — 55 ~^ 



